Collective behavior of squirmers in thin films.

Bacteria in biofilms form complex structures and can collectively migrate within mobile aggregates, which is referred to as swarming. This behavior is influenced by a combination of various factors, including morphological characteristics and propulsive forces of swimmers, their volume fraction within a confined environment, and hydrodynamic and steric interactions between them. In our study, we employ the squirmer model for microswimmers and the dissipative particle dynamics method for fluid modeling to investigate the collective motion of swimmers in thin films. The film thickness permits a free orientation of non-spherical squirmers, but constraints them to form a two-layered structure at maximum. Structural and dynamic properties of squirmer suspensions confined within the slit are analyzed for different volume fractions of swimmers, motility types (e.g., pusher, neutral squirmer, puller), and the presence of a rotlet dipolar flow field, which mimics the counter-rotating flow generated by flagellated bacteria. Different states are characterized, including a gas-like phase, swarming, and motility-induced phase separation, as a function of increasing volume fraction. Our study highlights the importance of an anisotropic swimmer shape, hydrodynamic interactions between squirmers, and their interaction with the walls for the emergence of different collective behaviors. Interestingly, the formation of collective structures may not be symmetric with respect to the two walls. Furthermore, the presence of a rotlet dipole significantly mitigates differences in the collective behavior between various swimmer types. These results contribute to a better understanding of the formation of bacterial biofilms and the emergence of collective states in confined active matter.

Motion of microswimmers in cylindrical microchannels.

Biological and artificial microswimmers often have to propel through a variety of environments, ranging from heterogeneous suspending media to strong geometrical confinement. Under confinement, local flow fields generated by microswimmers, and steric and hydrodynamic interactions with their environment determine the locomotion. We propose a squirmer-like model to describe the motion of microswimmers in cylindrical microchannels, where propulsion is generated by a fixed surface slip velocity. The model is studied using an approximate analytical solution for cylindrical swimmer shapes, and by numerical hydrodynamics simulations for spherical and spheroidal shapes. For the numerical simulations, we employ the dissipative particle dynamics method for modelling fluid flow. Both the analytical model and simulations show that the propulsion force increases with increasing confinement. However, the swimming velocity under confinement remains lower than the swimmer speed without confinement for all investigated conditions. In simulations, different swimming modes (i.e. pusher, neutral, puller) are investigated, and found to play a significant role in the generation of propulsion force when a swimmer approaches a dead end of a capillary tube. Propulsion generation in confined systems is local, such that the generated flow field generally vanishes beyond the characteristic size of the swimmer. These results contribute to a better understanding of microswimmer force generation and propulsion under strong confinement, including the motion in porous media and in narrow channels.

Dynamic shapes of floppy vesicles enclosing active Brownian particles with membrane adhesion.

Recent advances in micro- and nano-technologies allow the construction of complex active systems from biological and synthetic materials. An interesting example is active vesicles, which consist of a membrane enclosing self-propelled particles, and exhibit several features resembling biological cells. We investigate numerically the behavior of active vesicles, where the enclosed self-propelled particles can adhere to the membrane. A vesicle is represented by a dynamically triangulated membrane, while the adhesive active particles are modelled as active Brownian particles (ABPs) that interact with the membrane via the Lennard-Jones potential. Phase diagrams of dynamic vesicle shapes as a function of ABP activity and particle volume fraction inside the vesicle are constructed for different strengths of adhesive interactions. At low ABP activity, adhesive interactions dominate over the propulsion forces, such that the vesicle attains near static configurations, with protrusions of membrane-wrapped ABPs having ring-like and sheet-like structures. At moderate particle densities and strong enough activities, active vesicles show dynamic highly-branched tethers filled with string-like arrangements of ABPs, which do not occur in the absence of particle adhesion to the membrane. At large volume fractions of ABPs, vesicles fluctuate for moderate particle activities, and elongate and finally split into two vesicles for large ABP propulsion strengths. We also analyze membrane tension, active fluctuations, and ABP characteristics (e.g., mobility, clustering), and compare them to the case of active vesicles with non-adhesive ABPs. The adhesion of ABPs to the membrane significantly alters the behavior of active vesicles, and provides an additional parameter for controlling their behavior.

Competition between deformation and free volume quantified by 3D image analysis of red blood cell.

Cells in living organisms are subjected to mechanical strains caused by external forces like overcrowding, resulting in strong deformations that affect cell function. We study the interplay between deformation and crowding of red blood cells (RBCs) in dispersions of nonabsorbing rod-like viruses. We identify a sequence of configurational transitions of RBC doublets, including configurations that can only be induced by long-ranged attraction: highly fluctuating T-shaped and face-to-face configurations at low, and doublets approaching a complete spherical configuration at high, rod concentrations. Complementary simulations are used to explore different energy contributions to deformation as well as the stability of RBC doublet configurations. Our advanced analysis of 3D reconstructed confocal images of RBC doublets quantifies the depletion interaction and the resulting deformation energy. Thus, we introduce a noninvasive, high-throughput platform that is generally applicable to investigate the mechanical response of biological cells to external forces and characterize their mechanical properties.

Hydrodynamic clustering of two finite-length flagellated swimmers in viscoelastic fluids.

Clustering of flagellated microswimmers such as sperm is often mediated by hydrodynamic interactions between them. To better understand the interaction of microswimmers in viscoelastic fluids, we perform two-dimensional simulations of two swimming sheets, using a viscoelastic version of the smoothed dissipative particle dynamics method that implements the Oldroyd-B fluid model. Elasticity of sheets (stiff versus soft) defines two qualitatively different regimes of clustering, where stiff sheets exhibit a much more robust clustering than soft sheets. A formed doublet of soft sheets generally swims faster than a single swimmer, while a pair of two stiff sheets normally shows no speed enhancement after clustering. A pair of two identical swimmers is stable for most conditions, while differences in the beating amplitudes and/or frequencies between the two sheets can destroy the doublet stability. Clustering of two distinct swimmers is most stable at Deborah numbers of De = τω ≈ 1 (τ is the relaxation time of a viscoelastic fluid and ω is the beating frequency), in agreement with experimental observations. Therefore, the clustering of two swimmers depends non-monotonically on De. Our results suggest that the cluster stability is likely a dominant factor which determines the cluster size of collectively moving flagellated swimmers.

Correction: Non-equilibrium shapes and dynamics of active vesicles.

Correction for 'Non-equilibrium shapes and dynamics of active vesicles' by Priyanka Iyer et al., Soft Matter, 2022, 18, 6868-6881, https://doi.org/10.1039/D2SM00622G.

Non-equilibrium shapes and dynamics of active vesicles.

Active vesicles, constructed through the confinement of self-propelled particles (SPPs) inside a lipid membrane shell, exhibit a large variety of non-equilibrium shapes, ranging from the formation of local tethers and dendritic conformations, to prolate and bola-like structures. To better understand the behavior of active vesicles, we perform simulations of membranes modelled as dynamically triangulated surfaces enclosing active Brownian particles. A systematic analysis of membrane deformations and SPP clustering, as a function of SPP activity and volume fraction inside the vesicle is carried out. Distributions of membrane local curvature, and the clustering and mobility of SPPs obtained from simulations of active vesicles are analysed. There exists a feedback mechanism between the enhancement of membrane curvature, the formation of clusters of active particles, and local or global changes in vesicle shape. The emergence of active tension due to the activity of SPPs can well be captured by the Young-Laplace equation. Furthermore, a simple numerical method for tether detection is presented and used to determine correlations between the number of tethers, their length, and local curvature. We also provide several geometrical arguments to explain different tether characteristics for various conditions. These results contribute to the future development of steerable active vesicles or soft micro-robots whose behaviour can be controlled and used for potential applications.

Erythrocyte sedimentation: Effect of aggregation energy on gel structure during collapse.

The erythrocyte (or red blood cell) sedimentation rate (ESR) is commonly interpreted as a measure of cell aggregation and as a biomarker of inflammation. It is well known that an increase of fibrinogen concentration, an aggregation-inducing protein for erythrocytes, leads to an increase of the sedimentation rate of erythrocytes, which is generally explained through the formation and faster settling of large disjoint aggregates. However, many aspects of erythrocyte sedimentation conform well with the collapse of a particle gel rather than with the sedimentation of disjoint aggregates. Using experiments and cell-level numerical simulations, we systematically investigate the dependence of ESR on fibrinogen concentration and its relation to the microstructure of the gel-like erythrocyte suspension. We show that for physiological aggregation interactions, an increase in the attraction strength between cells results in a cell network with larger void spaces. This geometrical change in the network structure occurs due to anisotropic shape and deformability of erythrocytes and leads to an increased gel permeability and faster sedimentation. Our results provide a comprehensive relation between the ESR and the cell-level structure of erythrocyte suspensions and support the gel hypothesis in the interpretation of blood sedimentation.

Erythrocyte Sedimentation: Collapse of a High-Volume-Fraction Soft-Particle Gel.

The erythrocyte sedimentation rate is one of the oldest medical diagnostic methods whose physical mechanisms remain debatable today. Using both light microscopy and mesoscale cell-level simulations, we show that erythrocytes form a soft-particle gel. Furthermore, the high volume fraction of erythrocytes, their deformability, and weak attraction lead to unusual properties of this gel. A theoretical model for the gravitational collapse is developed, whose predictions are in agreement with detailed macroscopic measurements of the interface velocity.

Effect of malaria parasite shape on its alignment at erythrocyte membrane.

During the blood stage of malaria pathogenesis, parasites invade healthy red blood cells (RBC) to multiply inside the host and evade the immune response. When attached to RBC, the parasite first has to align its apex with the membrane for a successful invasion. Since the parasite's apex sits at the pointed end of an oval (egg-like) shape with a large local curvature, apical alignment is in general an energetically unfavorable process. Previously, using coarse-grained mesoscopic simulations, we have shown that optimal alignment time is achieved due to RBC membrane deformation and the stochastic nature of bond-based interactions between the parasite and RBC membrane (Hillringhaus et al., 2020). Here, we demonstrate that the parasite's shape has a prominent effect on the alignment process. The alignment times of spherical parasites for intermediate and large bond off-rates (or weak membrane-parasite interactions) are found to be close to those of an egg-like shape. However, for small bond off-rates (or strong adhesion and large membrane deformations), the alignment time for a spherical shape increases drastically. Parasite shapes with large aspect ratios such as oblate and long prolate ellipsoids are found to exhibit very long alignment times in comparison to the egg-like shape. At a stiffened RBC, a spherical parasite aligns faster than any other investigated shape. This study shows that the original egg-like shape performs not worse for parasite alignment than other considered shapes but is more robust with respect to different adhesion interactions and RBC membrane rigidities.

The Erythrocyte Sedimentation Rate and Its Relation to Cell Shape and Rigidity of Red Blood Cells from Chorea-Acanthocytosis Patients in an Off-Label Treatment with Dasatinib.

BACKGROUND: Chorea-acanthocytosis (ChAc) is a rare hereditary neurodegenerative disease with deformed red blood cells (RBCs), so-called acanthocytes, as a typical marker of the disease. Erythrocyte sedimentation rate (ESR) was recently proposed as a diagnostic biomarker. To date, there is no treatment option for affected patients, but promising therapy candidates, such as dasatinib, a Lyn-kinase inhibitor, have been identified. METHODS: RBCs of two ChAc patients during and after dasatinib treatment were characterized by the ESR, clinical hematology parameters and the 3D shape classification in stasis based on an artificial neural network. Furthermore, mathematical modeling was performed to understand the contribution of cell morphology and cell rigidity to the ESR. Microfluidic measurements were used to compare the RBC rigidity between ChAc patients and healthy controls. RESULTS: The mechano-morphological characterization of RBCs from two ChAc patients in an off-label treatment with dasatinib revealed differences in the ESR and the acanthocyte count during and after the treatment period, which could not directly be related to each other. Clinical hematology parameters were in the normal range. Mathematical modeling indicated that RBC rigidity is more important for delayed ESR than cell shape. Microfluidic experiments confirmed a higher rigidity in the normocytes of ChAc patients compared to healthy controls. CONCLUSIONS: The results increase our understanding of the role of acanthocytes and their associated properties in the ESR, but the data are too sparse to answer the question of whether the ESR is a suitable biomarker for treatment success, whereas a correlation between hematological and neuronal phenotype is still subject to verification.

Acanthocyte Sedimentation Rate as a Diagnostic Biomarker for Neuroacanthocytosis Syndromes: Experimental Evidence and Physical Justification.

(1) Background: Chorea-acanthocytosis and McLeod syndrome are the core diseases among the group of rare neurodegenerative disorders called neuroacanthocytosis syndromes (NASs). NAS patients have a variable number of irregularly spiky erythrocytes, so-called acanthocytes. Their detection is a crucial but error-prone parameter in the diagnosis of NASs, often leading to misdiagnoses. (2) Methods: We measured the standard Westergren erythrocyte sedimentation rate (ESR) of various blood samples from NAS patients and healthy controls. Furthermore, we manipulated the ESR by swapping the erythrocytes and plasma of different individuals, as well as replacing plasma with dextran. These measurements were complemented by clinical laboratory data and single-cell adhesion force measurements. Additionally, we followed theoretical modeling approaches. (3) Results: We show that the acanthocyte sedimentation rate (ASR) with a two-hour read-out is significantly prolonged in chorea-acanthocytosis and McLeod syndrome without overlap compared to the ESR of the controls. Mechanistically, through modern colloidal physics, we show that acanthocyte aggregation and plasma fibrinogen levels slow down the sedimentation. Moreover, the inverse of ASR correlates with the number of acanthocytes (R2=0.61, p=0.004). (4) Conclusions: The ASR/ESR is a clear, robust and easily obtainable diagnostic marker. Independently of NASs, we also regard this study as a hallmark of the physical view of erythrocyte sedimentation by describing anticoagulated blood in stasis as a percolating gel, allowing the application of colloidal physics theory.

Effect of cytosol viscosity on the flow behavior of red blood cell suspensions in microvessels.

OBJECTIVE: The flow behavior of blood is strongly affected by red blood cell (RBC) properties, such as the viscosity ratio C between cytosol and suspending medium, which can significantly be altered in several pathologies (e.g. sickle-cell disease, malaria). The main objective of this study is to understand the effect of C on macroscopic blood flow properties such as flow resistance in microvessels, and to link it to the deformation and dynamics of single RBCs. METHODS: We employ mesoscopic hydrodynamic simulations to investigate flow properties of RBC suspensions with different cytosol viscosities for various flow conditions in cylindrical microchannels. RESULTS: Starting from a dispersed cell configuration which approximates RBC dispersion at vessel bifurcations in the microvasculature, we find that the flow convergence and development of RBC-free layer (RBC-FL) depend only weakly on C, and require a convergence length in the range of 25D-50D, where D is channel diameter. In vessels with D≤20µm , the final resistance of developed flow is nearly the same for C = 5 and C = 1, while for D=40µm , the flow resistance for C = 5 is about 10% larger than for C = 1. The similarities and differences in flow resistance can be explained by viscosity-dependent RBC-FL thicknesses, which are associated with the viscosity-dependent dynamics of single RBCs. CONCLUSIONS: The weak effect on the flow resistance and RBC-FL explains why RBCs can contain a high concentration of hemoglobin for efficient oxygen delivery, without a pronounced increase in the flow resistance. Furthermore, our results suggest that significant alterations in microvascular flow in various pathologies are likely not due to mere changes in cytosolic viscosity.

Active particles induce large shape deformations in giant lipid vesicles.

Biological cells generate intricate structures by sculpting their membrane from within to actively sense and respond to external stimuli or to explore their environment1-4. Several pathogenic bacteria also provide examples of how localized forces strongly deform cell membranes from inside, leading to the invasion of neighbouring healthy mammalian cells5. Giant unilamellar vesicles have been successfully used as a minimal model system with which to mimic biological cells6-11, but the realization of a minimal system with localized active internal forces that can strongly deform lipid membranes from within and lead to dramatic shape changes remains challenging. Here we present a combined experimental and simulation study that demonstrates how self-propelled particles enclosed in giant unilamellar vesicles can induce a plethora of non-equilibrium shapes and active membrane fluctuations. Using confocal microscopy, in the experiments we explore the membrane response to local forces exerted by self-phoretic Janus microswimmers. To quantify dynamic membrane changes, we perform Langevin dynamics simulations of active Brownian particles enclosed in thin membrane shells modelled by dynamically triangulated surfaces. The most pronounced shape changes are observed at low and moderate particle loadings, with the formation of tether-like protrusions and highly branched, dendritic structures, whereas at high volume fractions globally deformed vesicle shapes are observed. The resulting state diagram predicts the conditions under which local internal forces generate various membrane shapes. A controlled realization of such distorted vesicle morphologies could improve the design of artificial systems such as small-scale soft robots and synthetic cells.

Deterministic Lateral Displacement: Challenges and Perspectives.

The advent of microfluidics in the 1990s promised a revolution in multiple industries from healthcare to chemical processing. Deterministic lateral displacement (DLD) is a continuous-flow microfluidic particle separation method discovered in 2004 that has been applied successfully and widely to the separation of blood cells, yeast, spores, bacteria, viruses, DNA, droplets, and more. Deterministic lateral displacement is conceptually simple and can deliver consistent performance over a wide range of flow rates and particle concentrations. Despite wide use and in-depth study, DLD has not yet been fully elucidated or optimized, with different approaches to the same problem yielding varying results. We endeavor here to provide up-to-date expert opinion on the state-of-art and current fundamental, practical, and commercial challenges with DLD as well as describe experimental and modeling opportunities. Because these challenges and opportunities arise from constraints on hydrodynamics, fabrication, and operation at the micro- and nanoscale, we expect this Perspective to serve as a guide for the broader micro- and nanofluidic community to identify and to address open questions in the field.

Stochastic bond dynamics facilitates alignment of malaria parasite at erythrocyte membrane upon invasion.

Malaria parasites invade healthy red blood cells (RBCs) during the blood stage of the disease. Even though parasites initially adhere to RBCs with a random orientation, they need to align their apex toward the membrane in order to start the invasion process. Using hydrodynamic simulations of a RBC and parasite, where both interact through discrete stochastic bonds, we show that parasite alignment is governed by the combination of RBC membrane deformability and dynamics of adhesion bonds. The stochastic nature of bond-based interactions facilitates a diffusive-like re-orientation of the parasite at the RBC membrane, while RBC deformation aids in the establishment of apex-membrane contact through partial parasite wrapping by the membrane. This bond-based model for parasite adhesion quantitatively captures alignment times measured experimentally and demonstrates that alignment times increase drastically with increasing rigidity of the RBC membrane. Our results suggest that the alignment process is mediated simply by passive parasite adhesion.

Dissipative particle dynamics with energy conservation: Isoenergetic integration and transport properties.

Simulations of nano- to micro-meter scale fluidic systems under thermal gradients require consistent mesoscopic methods accounting for both hydrodynamic interactions and proper transport of energy. One such method is dissipative particle dynamics with energy conservation (DPDE), which has been used for various fluid systems with non-uniform temperature distributions. We propose an easily parallelizable modification of the velocity-Verlet algorithm based on local energy redistribution for each DPDE particle such that the total energy in a simulated system is conserved up to machine precision. Furthermore, transport properties of a DPDE fluid are analyzed in detail. In particular, an analytical approximation for the thermal conductivity coefficient is derived, which allows its a priori estimation for a given parameter set. Finally, we provide approximate expressions for the dimensionless Prandtl and Schmidt numbers, which characterize fluid transport properties and can be adjusted independently by a proper selection of model parameters. In conclusion, our results strengthen the DPDE method as a very robust approach for the investigation of mesoscopic systems with temperature inhomogeneities.

Importance of Erythrocyte Deformability for the Alignment of Malaria Parasite upon Invasion.

Invasion of erythrocytes by merozoites is an essential step for the survival and progression of malaria parasites. To invade red blood cells (RBCs), apicomplexan parasites have to adhere with their apex to the RBC membrane. This necessary apex-membrane contact (or alignment) is not immediately established because the orientation of a free merozoite with respect to the RBC membrane is random when an adhesion contact first occurs. Therefore, it has been suggested that after the initial adhesion, merozoites facilitate their proper alignment by inducing considerable membrane deformations, frequently observed before the invasion process. This proposition is based on a positive correlation between RBC membrane deformation and successful invasion; however, the role of RBC mechanics and its deformation in the alignment process remains elusive. Using a mechanically realistic model of a deformable RBC, we investigate numerically the importance of RBC deformability for merozoite alignment. Adhesion between the parasite and RBC membrane is modeled by an attractive potential that might be inhomogeneous, mimicking possible adhesion gradients at the surface of a parasite. Our results show that RBC membrane deformations are crucial for successful merozoite alignment and require interaction strengths comparable to adhesion forces measured experimentally. Adhesion gradients along the parasite body further improve its alignment. Finally, an increased membrane rigidity is found to result in poor merozoite alignment, which can be a possible reason for a reduction in the invasion susceptibility of RBCs in several blood diseases associated with membrane stiffening.